TPI-MINN-97/30-T 
NUC-MINN-97/14-T 
HEP-MINN-97-1614 
LBNL 41103 



Temperature Dependence of 77 and rf Masses 

o 

(N 

> ! Jamal Jalilian-Marian 1 and Bayram Tekin 2 

Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720 
2 School of Physics and Astronomy, University of Minnesota, Minneapolis, MN 55455 

r-> ' 

^3 
a: 

> : 

•'-i ■ Abstract 

X 

■ We investigate the temperature dependence of n and rf masses due to scatterings 

from thermal pions in a heat bath using the non-linear sigma model. We show that 
mass shifts of n and rf and the shift in the mixing angle are negligible. 



1 Introduction 



Temperature dependence of r\ and vj masses is a very interesting subject and has been 
a hot topic lately Jl|, ^ The rj', being heavier than the other pseudoscalar mesons, 
naturally begs the question why it is so heavy. Besides being of academic interest, this 
also has important experimental consequences: for instance, in dilepton spectra which may 
help with identifying signals of a possible quark-gluon plasma. 

Even though one can understand the origin of this mass difference between the 7/ 
and the other pseudoscalars at T = in terms of the U(l) anomaly and instantons, a 
satisfactory understanding of this problem at finite temperatures is still missing. On one 
hand, one might expect that at high enough temperatures instanton effects and the anomaly 
will go away and as a result rf will become mass degenerate with 77 and other pseudoscalar 
mesons [0. The first question that comes to mind is how high does the temperature have 
to be since the relavent temperature range of 100 — 200 MeV may not be high enough in 
order to ignore instanton effects. On the other hand, it has been argued that instantons 
do not go away in the temperature range of interest but rather are rearranged so that the 
effect of the U(l) anomaly is essentially unchanged. This leads to the non-strange meson 
being heavier than the strange one if the anomaly is strong enough Q. 

In this work, our goal is rather modest and pedagogical. We propose to use the non- 
linear sigma model, with symmetry breaking and an anomaly term added, to investigate 
the temperature dependence of 77 and rf masses for temperatures up to 150 MeV. We 
consider the change in masses of rj and rf due to scatterings from a heat bath consisting 
mainly of pions. By evaluating the relevant one loop diagrams at finite temperature, one 
can calculate the temperature dependence of this change in masses. 

One may question the validity of using the non-linear sigma model for temperatures 
around the chiral phase transition temperature T ~ 150 — 200. However, we are rather 
interested in gross features of the model with respect to the mass shifts below the phase 
transition temperature and certainly do not claim to address the problem of mass change 
across the phase transition temperature. As a first approximation, this will tell us the 
direction and magnitude of mass shifts. 
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Furthermore, we will neglect the effect of scatterings from other particles in the heat 
bath like kaons, etc. since they are much heavier than pions and one does not expect 
to have too many of them in the heat bath in temperature ranges we are interested in. 
We will take the pion mass to be constant. This is certainly a good approximation for 
temperatures up to 100 MeV but will start changing beyond that where scattering from 
other particles may also become important. This is all ignored. 

Another problem, and perhaps the central one is the temperature dependence of the 
coefficient of the anomaly ||, ^ 0, ||, |9| . The temperature dependence of this coefficient 
is not well known and is usually modeled in the literature such that it goes to zero at 
sufficiently high temperatures. Here, we will take this coefficient to be constant. This is 
plausible in the temperature region we are interested as can be seen, for instance, in Fig. 
(1) of the paper by Schafer || where it is shown that the t'Hooft operator varies only slowly 
up to T ~ 100 MeV and that only around T ~ 150 MeV does it have a sharp increase. 
Again, we emphasize that it is not our goal in this work to investigate the behavior of 
masses around the chiral phase transition temperature. 

We then consider the one loop diagrams with pions in the loop. These diagrams, 
evaluated at finite T, correspond to scattering of rj or rf from thermal pions in the heat 
bath. We show that the changes in masses are negligible. We then consider the possibility 
that 7] or rf can fluctuate into a vector or axial vector meson, scatter from thermal pions 
in the bath and then go back to its original state. We show that these processes are absent 
and do not contribute to mass shifts of rj and rf . 

2 Pseudoscalars 

We will start with the following effective Lagrangian describing the pseudoscalar nonet, 
£ = ^tid^Ud^ + If^tM(U + U ] -2) + a(detU + detU 1 -2), (1) 
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where U = exp (f^) and F w = V2f n = 132 MeV is the pion decay constant. The mass 



matrix for the nonet is 
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The first two terms in our effective Lagrangian are the usual non-linear sigma model with 
massive psuedoscalar fields given by 

f -1=^0 + 2^0) n + K + \ 

r _i_ Vs+V^rio 
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We have included a determinantal term in our effective Lagrangian to account for the 
axial U{\) anomaly. The coefficient of the anomaly term, a, is to be determined from 
experimental data at zero temperature. With this effective Lagrangian at hand, we will 
expand the field U and look for interaction terms between rj and rj' and pions. We ignore 
kaons altogether because they are much heavier than pions. In other words, our heat bath 
consists mainly of pions. 

The quadratic piece of the Lagrangian then becomes 



1 i 
£ quad = 2^ o)2 + o^ m ^ 
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whereas the interaction terms are 
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It is interesting to note that there are no derivative coupling involving 77 , r( with pions. 

At the classical (tree) level, it is clear from the quadratic part of our Lagrangian that 
we need to rotate the 770 and r] s fields and write them in terms of the physical fields i] and 
7/ defined as 

i] = —r]osm9 + r]$cos9, 



rjo cos 9 + r] 8 sin 9. 



(6) 



The mixing angle 9 is defined as 
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The physical masses of 77 and rj are given by 



m 2 , 



C 



B cos 9 + A sin sin 20 
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5 cos 2 9 + A sin 2 + — sin 2# 



We determine the mixing angle and the coefficient of the anomaly by matching the 77 
and 7/ masses with the experimental values. Using m n = 135 MeV, mx = 497 MeV and 
a = |j- [m^+m 2 — 2m|-], the mixing angle comes out to be # = —18.5° in close agreement 
with the experimental value of —20°. This redefinition would then diagonalize the quadratic 
part of the Lagrangian in the usual manner. However, here we will be considering quantum 
corrections and have to include the quantum loop effects coming from the pion tadpole 
before diagonalizing the mass matrix. 

In order to calculate the mass shifts due to scattering with thermal pions, all we need 
to do now is to construct the lowest order Feynman diagrams contributing to the self energy 
matrix Y\. The components of this matrix are represented diagrammatically by 





n = _ vs \ y vo 

Diagram 1: Relevant tadpoles for the self energy matrix. 
These diagrams give 

]J = -2A (1 , 2 , 3 )T£ / f P * (9) 

where 



Since we are interested in the temperature dependence of these corrections, we will 
neglect the vacuum part of the self energy diagram and assume that it is renormalized as 
in zero temperature field theory ||10|| . The T dependent pieces of the self energy diagram 
will then give the change in the masses of 770 and rj 8 . They are proportional to the integral 

d 3 p 1 1 



(2tt) 3 



(11) 



where (3 = 1/T and w = yp 2 + m 2 . This integral can be evaluated analytically in terms 
of an infinite sum of modified Bessel functions K\ (x) jlT| with the result 



1 = ^1^ A (12) 
2tt 2 T 

We are interested in not too high temperatures so that we will take m, = 135 MeV to 
be the upper limit for T. One can then approximate the exact series by the asymptotic 



limit jnj where, for x > 1, 
K v {x) = ^e~ x 



1 | (4^ 2 -l 2 ) (4^ 2 -l 2 )(4z, 2 -3 2 ) 
l!8x 2!(8a;) 2 



(13) 



One can check numerically that it is enough to keep the first three terms in the infinite 
series. It turns out to be a good approximation even when x ~ 1 as will be shown below. 
Keeping the first three terms, we have 



2?r V 2vr 



3T 15T 2 1 3T 15T 2 , 

1 H 1 e~~(l H ) 

8m, 128m 2 2^ 16m, 512m 2 ' 

2 



1 2m„ 3T 15T 

+ — =e"— 1 + 



3v/3 v 24m, 1152m 2 / 

6 



(14) 



To check the validity of this approximation near x = 1 where our approximation would 
have the largest error , one can evaluate the value of integral / exactly at this value. One 
obtains / = 0.035m^. Equation ( |14D differs from this value by 2 percent. 

Including these corrections in the original Lagrangian (|J) we get 
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Now we have to diagonalize the mass matrix in order to get the r\ and 7/ masses. Define 
the physical fields rj and 7/ in terms of r/o and r/g as 



77 = —770 sin + rjs cos < 



r/o cos + 773 sin 0. 



(16) 



The new mixing angle 0, defined as 



tan 20 



C 



A-B 



where the constants A, B, C are now 



A = A-^I 

2ml 

B = B fl 



(17) 



The mixing angle at finite temperature is not significantly different from its value at zero 
temperature. 

At finite temperature we obtain the following formulae for the masses: 

TYl? i 

m 2 v (T) = m 2 (o) __^ (1 + sin 2 ___ sin20)/ 

rri ^ j 

m 2 v ,(T) = m 2 ;(o) __^ (1 + cos 2 0+ __ sin20)L (18) 
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Figure 1: function of temperature. 
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Figure 2: m v > as a function of temperature. 
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It is interesting to note that the correction to r\ mass is larger than the correction to mass 
of rf . Since the integral I is very small even at 135 MeV it is clear that neither mass 
receives a significant change. In Figures 1 and 2, we plot the values of 77 and rj' mass as a 
function of temperature. 

One may naturally wonder if inclusion of vector and axial vector mesons would change 
our results drastically. For instance, interaction of the type depicted in Diagram (2) can 
contribute to the change in the 77, r( masses. 



710 




V,A 

Diagram 2: Vector/axial vector contribution. 

To investigate such a possibility, we will need to write an effective Lagrangian which in- 
cludes both vectors and axial vectors in addition to our pseudoscalar nonet. We will 
start with the following effective Lagrangian introduced by Meissner|T^[] based on massive 
Yang-Mills gauge theory 



£ = C + Ci + £ 2 + £3 

where 

A, = ltr[(D,U)(D,U)\ 

£2 = mlti[A L A^ L + A R A^ R ] + Bti[A L UA^ R U ] ], 

£ 3 = If^tM(U + U 1 -2) + a(detU + detU^ -2). (19) 
8 

The covariant derivative is defined as 

D^U = d,-igA L ^+igA R (20) 
9 



where the field strength is 



Ftf = d»A L * - d u A L ' R - ig[A L ' R , A L ^ R ] (21) 



and left-handed and right-handed fields A^' R are defined as 



Ai=\{Vv + Ap), A R = l -{V,-A,). (22) 



Notice that at this point both vector and axial vector mesons have the same mass m 2 , and 
7 and B are related to m p [p~2j] . 

In order to avoid the difficulties related to experimental uncertainties with regard to 
axial vector mesons, one can eliminate the explicit dependence of the Lagrangian on axial 
vectors by a suitable gauge transformations such that all axial vectors are zero in this new 
gauge as done, for example, in [O. We then have 



C x = ( 7 -l)tr[F^(p)F^(p)], 

C 2 = (B + 2m 2 )tr[p^] + (Vs)(£ + 2m 

+ 2(ml/g 2 )ti[d, l U 1/2 d^U- 1/2 } - (B / g 2 )ti[U- 1/2 d, t U- 1/2 U 1/2 d' l U 1/2 } (23) 

Choosing 7 = 3/4 will give the correct kinetic term for the vector mesons. Also, m 2 = 
4m 2 , + 2B. For a discussion of these parameters see [|12|]. In order to get the appropriate 
interaction terms for pseudoscalars and vector mesons we expand the Lagrangian in terms 
of the pseudoscalar fields as before and examine the relevant vertex, which is of the form 

trV> ^ = V^%<p c tr\ a [A 6 , A c ] (24) 

where A = 1 while for a = 1, • ■ ■ 8 they are the usual Gell-Mann matrices. Since we 
are interested in interactions of r] and t] 8 (corresponding to A and A 8 ) and pions (A 1 ' 2,3 ) 
with vectors, then a,b,c = 0, 1, 2, 3, 8 only. It is then easy to show that this term does 
not involve interaction vertices of the type rjirp which would contribute to the change in 
masses of rj and rj' through diagrams depicted in Diagram (2). 
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3 Conclusion 



We used the non-linear sigma model, with appropriately added symmetry breaking and 
anomaly terms, to investigate the shifts in 77 and rj masses due to scatterings with thermal 
pions in a heat bath. We showed that the mass shifts as well as the change in the mixing 
angle are negligible. We also considered possible effects of vector and axial vector mesons 
on 77 and 7/ masses and showed that vector and axial vector mesons do not contribute. 
Our straight-forward calcualtion imply that if experiments show a stronger temperature 
dependence , then one would conclude that non-linear sigma model description fails to be 
true and one should look for new physics. 
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